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the bilayer by 10% . 
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1. INTRODUCTION 

Heterogeneous ferromagnetic-superconducting (FM- 
SC) systems have attracted much attention recently 0. 
In majority of these systems the proximity effect is sup- 
pressed by the oxide layer between FM and SC. Inhomo- 
geneous magnetization produces a magnetic field pene- 
trating into the superconductor and inducing supercur- 
rents. Supercurrents in turn produce the magnetic field 
acting on the magnetization. Thus, the FM and SC of 
systems strongly interact via magnetic field. Systems in 
which both, FM and SC parts are thin films represent a 
special interest for the experiment and can be analyzed 
theoretically. In these systems, spontaneous vortices ap- 
pear due to the magnetic interaction 0. Erdin et al. Q 
have developed a method to calculate the arrangement 
of the magnetization in the FM film and supercurrents 
including vortices in the SC film in the London's approx- 
imation. The London's approximation is quite good for 
these mesoscopic systems because characteristic length 
scales for magnetic field (the effective penetration depth, 
the linear size of the textures) are much larger than the 
coherence length £ of the superconductor. This method 
was applied recently Q to study topological textures in 
a FM-SC bilayer (FSB). It was shown that the homo- 
geneous state of the FSB with the magnetization per- 
pendicular to the layer is unstable with respect to the 
formation of vortices. The ground state of the FSB rep- 
resents a periodic array of stripe domains in which the 
direction of the magnetization in the FM film and the 
vorticity in the SC film alternate together. 

In this article we study how the combined FM-SC tex- 
tures change the superconducting transition temperature 
in bilayers and multilayers. For this purpose we extend 
theory of spontaneous SC-FM structures developed in 
|4| to the case of multilayers. We find that spontaneous 
domain- vortex structures increase the transition temper- 
ature, whereas vortex structures induced by a periodic 
array of magnetic dots decreases the transition tempera- 
ture. The magnitude of this effect is large enough to en- 



able its experimental observation. Though the influence 
of the textures on the transition temperature is akin to 
the influence of the homogeneous magnetic field, there 
are important differences between these two phenomena: 
first, the average magnetic field may be zero for mag- 
netic textures; second, the reciprocal action of the mag- 
netic field generated by vortices onto magnetization is 
substantial. 

The plan of this article is as follows. In the next sec- 
tion we consider the change of the transition temperature 
in the spontaneous stripe structure formed in the FSB. 
In Section 3 we analyze how this stripe structure and 
the transition temperature change in the presence of an 
external magnetic field. In Section 4 we study the shift 
of the transition temperature in a square array of mag- 
netic dots. Section 5 is devoted to theory of the FM-SC 
spontaneous textures in multilayers and to the transition 
temperature shift in them. Our conclusion are given in 
Section 6. 



2. TRANSITION TEMPERATURE IN 
SPONTANEOUS STRIPE STRUCTURE OF FSB 

As it was mentioned earlier, the homogeneous state 
of the FSB with the magnetization perpendicular to the 
layer is unstable with respect to formation of a stripe 
domain structure, in which both, the the direction of 
the magnetization in the FM film and the circulation of 
the vortices in the SC film alternate together. Let the 
stripe width be L s . The magnetization can be written 
as m = ms{x)z, where the coordinate x is along the 
direction perpendicular to the domain walls, z denotes 
the unit vector perpendicular to the layers, s{x) is the 
periodic step function with period 2L S : 




+ 1 < x < L s , 
-1 L s < x < 2L S . 



The energy of the stripe structure per unit area U and 
the stripe equilibrium width L s were calculated in 0]. 
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Here we correct a calculational mistake of this work || 
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The notations in eqs. QJ, (J5J) are as follows: A e = ^ 
is the effective penetration depth in the SC him, whose 
thickness is denoted d s ; e dw is the linear tension of the 
domain wall; m — m — e v /(f>o] m is the magnetization 

■ In ^ is the 



per unit area of the FM film; e„ = 16 ^° x ^ — 
single vortex energy in the absence of the FM film; C — 
0.57721 • • • is the Euler constant. To find the transition 
temperature, we combine the energy given by eq. l(T)l 
with the Ginzburg-Landau free energy. The total free 
energy per unit area reads: 



F = U + F GL = 



-16m 2 -e d 
— exp( 



4m 2 



•C-l) 
-n s d s [a(T - T c ) + ^n s ] . 
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Here a and /3 are the Ginzburg-Landau parameters. We 
omit the gradient term in the Ginzburg-Landau equa- 
tion since the gradients of the phase are included in the 
energy (JIJ, whereas the gradients of the superconducting 
electrons density can be neglected everywhere beyond the 
vortex cores. Recalling that X\ = 4 ™^ c e a and plugging it 
into eq. QJ, we find: 
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where xo 
over n 
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and using the condition that n. 
transition temperature T* we obtain: 



Minimizing the total free energy 
at new 



AT r = T* - T r = 
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+ C-1). (5) 



Eq. (0 demonstrates that the interaction between FM 
and SC layers in the spontaneous stripe structure in- 
creases the transition temperature T c . Theory |4j as- 
sumes that the ratio tdw/m 2 is larger than 1, so that 
exp(|4^) ^> 1. However, this ratio can not be too large. 
Otherwise the width of domain becomes larger than the 
sample linear size. The maximum possible shift of tran- 
sition temperature corresponds to ~ 1. It is 
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The value a can be estimated as a ~ T c /ep, where 
cf is the Fermi energy. It is about 1CP 4 for low- 
temperature superconductors and about 1CP 3 — 10~ 2 for 
high-temperature superconductors. If we take AttM ~ 
IT, T c ~ 3K, e F ~ 30, 000K and d m ~ 300A, we obtain 
AT C /T C ~ 0.1. The dependence AT C oc m 2 = M 2 d 2 m on 
the thickness of the FM film and magnetization can be 
checked experimentally. 



3. SPONTANEOUS STRIPE STRUCTURE IN 
EXTERNAL FIELD 

In this section we study the spontaneous stripe system 
in the FM-SC bilayer in the presence of an external per- 
pendicular magnetic field B (along z direction). Since 
the magnetization tends to align along the external mag- 
netic field we anticipate that the width L\ of stripes with 
the magnetization parallel to the external magnetic field 
increases, whereas the width Li of the stripes with the 
antiparallel magnetization decreases. Let us define a step 
function with the period L = L\ + L% as follows: 



s{x) 



-1 (0<x<Li), 
-1 (£i < x < 2L 2 ) 



The Fourier transform of s(x) is: 



sg 



2i(l 
(ii- 



JGL, 



L 2 )/L 



)/{LG) (G^O), 
(G = 0). 



(7) 



Here G = 2irr/L and r are integers. For the sake of 
brevity, we denote t = L\ — L%, At large distance from the 
bilayer the magnetic field asymptotically becomes equal 
to the external magnetic field. The total magnetic flux 
is the same in any cross- section of the space. Thus, the 
average magnetic field through superconducting layer is 



— [ n(x)(f)odx = B e 
L Jo 



(8) 



This constraint can be incorporated by the standard La- 
grange multiplier method. According to it, the prob- 
lem is reduced to minimization of the effective energy 
U v — U v + U u , where U v is the energy of the vortex 
system, and U u = lo^- J n(x)x, to is the Lagrange multi- 
plier. The general expression for the free energy of a peri- 
odic stripe system of magnetization and vortices is given 
by equation (10) of the work 0]. Employing this equation 
and the Fourier expansion for the step function s(x) (see 
equation (JJJ) and denoting ng the Fourier-transform of 
the vortex density n(x), we obtain: 

Uv = ^s G n-G + 5 ^2 V G n G n^ G 
G G#o 
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where V G — ^/(2tt|G|) is the Fourier-transform of the 
vortex interaction. An infinitely large interaction term 
Vq—qTIq—qTIq—q has been omitted since it corresponds to 
the fixed average magnetic field. From equation ( |5J we 
readily find that the constraint condition implies: 



B„ 



nc=a 



(10) 



This equation confirms that Vg=o"G=o^G=o is the energy 
of the uniform external field. Minimization of the total 
vortex energy U v over the vortex density nc results in a 
system of equations: 



e v s G + V G n G = (G^O), 
e v s G = = Afc (G = 0). 



Plugging the solutions of ug(G ^ 0), u> and n G =o from 
equation (|1U|> . we finally obtain 

U = ^[G + m~ + ~ln(2 + 2cos^)] 



o 
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We are now in a position to minimize the total energy U 
over L and t. After doing that, we obtain two equations: 

C - 1 + In ^ - ln(l + tan 2 u) 1 / 2 - % = ,(12) 



tanu 



47TTO 



(13) 



Here u = £t. From equations (|12|l and l|13[) we find 
solutions for L and t: 



L = 



2U 



1 _ ( L s B c!ct \2 
± \ 2irm > 



2L LB ext 

t = — arctan — 
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(14) 



(15) 



where L s is given by equation (J2J. The results. i|T4*|) 
and (|15l) are similar to those for a purely ferromagnetic 
stripe structure in a single ferromagnetic film |l0|. The 
critical external field B c ext at which the domain structure 
vanishes is 



b: 



2irm/L s , 



(16) 



which is in the range of 1 — 10 Oe. 



In conclusion of this section, we consider how the 
temperature of superconducting transition of the bilayer 
changes in the presence of external magnetic field. Since 
at the field B° xt ~ 1 — 10 Oe the stripe structure vanishes, 
the superconducting transition proceeds in the homoge- 
neous state of ferromagnetic film excluding very small 
vicinity of T c . Therefore, it is determined by the nucle- 
ation process as in the case of a single superconducting 



film . The nucleation in a thin film for the field peren- 
dicular to it was considered by Tinkham [l^. Though 
the geometry is different from the bulk geometry consid- 
ered by Abrikosov^^, his solution can be directly ap- 
plied. The order parameter coincides with the Landau 
wave function for the first Landau level. In the case of 
bilayer the energy of the nucleus reads: 



U 



1 

2m 



-V A 

i c 



d 2 x 



AU. 



It differs from the energy in the absence of magnetic film 



by the value AU = —m J Bz'cPx, where Bj 1 is the 
magnetic field generated by the nucleus at the ferromag- 
netic film. The magnetic field generated by the nucleus 
reads: 



B(»)(i 



1 



- x J„(r')dV . (18) 



We assume that the current flows in the x — y plane. 
Since it has zero divergence, it can be represented as 

j„ =z xV/ , where f(x,y) is a function localized in a 
finite part of the plane. The flux of the induced field is: 



B^d 2 x 
1 



Z XV: 



z xV/(r') (fxcfx 1 (19) 



Simple transformations turn this integral into a following 
form: 



B ( ^d 2 x = ~\\ /(r')V 2 * r d 2 xd 3 x' . (20) 



This integral is equal to zero. Indeed, the 3d Laplasian 
of the inverse distance is proportional to 5 (r — r'). The 
radius-vector r belongs to the magnetic layer, whereas 
the radius-vector r' belongs to the superconducting one. 
What stays in equation (|20J) is the 2d Laplasian operator, 
which differs from the 3d one by the second derivative 
(~§z) \r-r'\ ' l ne integral from this function over the 
plane r is equal to zero when r' does not belong to the 
same plane. 

Thus, the interaction of the superconducting nucleus 
and the homogeneously magnetized film is zero indepen- 
dently on the wave function of the localized nucleus. 
Therefore, the transition temperature is the same as in 
the absence of the ferromagnetic film. Earlier we have 
found that, without magnetic field the ferromagnetic 
layer increases the transition temperature and its change 
may be about 1/10 of initial value. Surprisingly this 
shift disappears at extremely small external field about 
1-10 Gs. The reason of this paradox consist of unlimited 
expansion of domains upon approaching the transition 
point. An interesting and counterintuitive result of this 
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consideration is that a very small magnetic field about 
lOOe is enough to change the transition temperature rel- 
atively by 1/10. After this fast change the transition 
temperature changes substantially when magnetic field 
becomes comparable to H C 2, normally in the range IT. 



4. TRANSITION TEMPERATURE IN A 
SUPERCONDUCTING FILM WITH A SQUARE 
ARRAY OF FERROMAGNETIC DOTS 

Recently Erdin considered theoretically the vortex- 
antivortex textures in a superconducting film supplied 
with a regular square lattice of ferromagnetic dots [llj . 
Unfortunately, there occurred a mistake in his code for 
numerical calculations. Therefore, we reproduce here a 
part of his analysis, but our conclusions are completely 
different. We denote the dot lattice constant be a and 
assume each dot to be a circular thin disk with the radius 
R and constant magnetization m per unit area directed 
perpendicular to the plane (along z-axis). The total en- 















FIG. 1: Schematics representaion of FM dots with sponta- 
neous vortices and antivortices. The circles drawn by solid 
line represent FM dots. The dash half-circles with clockwise 
and anticlockwise arrow indicate vortices and antivortices re- 
spectively. 



ergy per unit area of the system is (ll| |: 

U Uijv T Umu T Um.m. • 



(21) 



The three terms in the right-hand side of the above equa- 
tion can be written as follows: 



A-kcl 1 G(l + 2AG) 

- mzGF ~ G 
a 2 ^ 1 + 2AG ' 

= ~ 2 - A E 1 + 2AG • 

G 



(22) 
(23) 
(24) 



where G = ^(r, s) (r, s are integers) are the reciprocal 
lattice vectors; Fg — J2i nie lG ' ri is the structure factor 
of the vortex lattice; n,i, indicate the vorticity and the 
position of a i-th vortex. Only the change of energy u mm 
in superconducting state compared to the normal state 
matters: 

^mm — Umm{^) W mm (A ► Oo) (^5) 

The last term in the r.-h.s. of equation (|25|l is the dipolar 
energy of the FM dots above the superconducting transi- 
tion. Below the superconducting transition temperature 
T c the magnetic field generated by the dots penetrates 
into the SC film and can create vortices and antivortices 
if the magnetization and the size of the dots are large 
enough [J- Erdin ^l| considered the case when only one 
vortex and one antivortex appear per a magnetic dot. He 
stated that there is a symmetry violation in the lowest 
energy state. Our numerical results contradict to this 
statement: we find that the vortex centers are located 
precisely under the centers of the magnetic dots, whereas 
the antivortex centers are located at the centers of the 
magnetic dot lattice. Employing this fact and keeping in 
mind that A > L near the new transition temperature 
T*, we can rewrite the total energy (|2*T|l as follows: 



4> e 2 d s n s a 
27tm s c 2 a 2 £ 



^d 2 n 2 s 
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where ^ means that the term r 
ii, I2 and -Z3 are defined as series: 
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are omitted. 
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(27) 



We combine this energy with the Ginzburg-Landau free 
energy for the SC film as it was done in the stripe struc- 
ture case: 



F = 



<p e 2 d s n s a (/)f,e 4 d 2 s n 2 s 
2irm s c 2 a 2 £ .>■<■■ /: ' 



t>le 2 d s n s 
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The condition of minimum over n s fro the free energy 1251 
reads: 

4> e 2 d s a §\e^d 2 s n s 

In - — -^—^ — s— <— io 



2i:m s c 2 a 2 £ 2ir 2 m 2 c 4 a 

<ble 2 d s , 4ir 2 mR s 2i: 2 m 2 e 2 d s R 2 
ih + — ; h) + ^—h 



Air 2 a 2 m s c 2 (po 
+a(T — T c )d s + /3n s d s = 



(29) 



At a new critical temperature T* the density of supercon- 
ducting carriers must be zero. Plugging n s (T*) = into 
eq. I|29l) . we obtain the shift of the critical temperature: 

AT r = . J° — ( ; h+h 



A-K 2 m s c 2 a 2 (pa 

a 8ir A m 2 R 2 
-27rln- -J 



h) ■ (30) 



Figure (2) shows the relation between AT C and R respec- 
tively for £ = 0.21a. To ensure spontaneous occurrence 
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FIG. 2: AT C vs. R for £ = 0.21a respectively for r = 10.0, 
12, 5 and 15.0, here r = 47r ' mL . AT C is in the unit 



of the vortices the inequality u mlJ + u vv < must be 
satisfied. It is equivalent to the following relation: 

4:TT 2 mR . 



-h +h - 2?rln- < 0. 



(31) 



The validity of the London's approximation implies £ <C 
a. This condition is violated in a close vicinity of the 
transition temperature, the smaller the larger is the dot 
lattice constant a. For a ~ lpm this vicinity is less or 
of the order of 0.01T C and further we neglect it. Fig. 
(2) shows that the shift of transition temperature is a 
rather complicated function of the dots radius R and the 
ratio r = 4ir 2 mL/(f>o. For each value r, there exists a 
threshold radius Ro, at which the vortices first appear. 
The shift of the transition temperature grows near the 
threshold with R/a until maximum and then decreases, 
it remain negative in the interval Ro/a and 1/2. For 
each fixed R/a, the shift of the transition temperature 
increases with the ratio r and is negtive. 



5. FERROMAGNETIC TEXTURES IN 
MULTILAYERS 

Let us with a FM-SC multilayer consisting of N bilay- 
ers, each having the thickness d in the limit Nd L s , 
where L s is the lateral size of the layers. If the mag- 
netic films are magnetized perpendicularly, the average 
induction inside the multilayer is B = Amvnj d and it is di- 
rected perpendicularly to the layers. The situation is the 
same as in the layered superconductors placed into an ex- 
ternal magnetic field 0. Therefore, pancake vortices in 
each superconducting layer must appear. Together they 
form the Abrikosov linear vortices and satisfy a condition: 
m<pQ jd > e; , which guarantees that the vortex line is en- 
ergy favorable. Here e; = e m j 01 is the usual vortex 
line energy per unit length, eo = ^o/(4 7r ^) 2 an d A is the 
planar bulk penetration depth. There is no need to con- 
sider the Josephson coupling effect in this case since the 
phase difference between superconducting layers is zero 
if the vortex lines are perpendicular to the layers. On the 
other hand, the Josephson vortices appear along the lay- 
ers if the magnetization m is parallel to the layers and 
satisfy a condition mcpo/d > ej where ej — 760 In is 
the Josephson vortex line energy and 7 is the anisotropy 
parameter for the layered superconductor 0. These ideas 
were applied by M. Houzet et al. Q to explain the mag- 
netic properties of the magnetic superconductor RuSr2 
GdCu 2 Os. In this article we presumably study the op- 
posite limit Nd <C A, where A = X 2 /d is the effective 
penetration depth for layered superconductors. 

In this section we first focus on few-layer superconduc- 
tors without ferromagnetic texture, which will be dis- 
cussed later. Pancake vortices in a finite stack of layers 
were discussed by Mints et al. 01 • We modify a method 
developed by Efetov for a plane superconductor and 
by K. Fischer ^| for a layered superconductor with in- 
finite number of layers. We consider a superconductor 
consisting of N layers coupled only by their magnetic 
field. To simplify the calculation, we assume that layers 
are infinitely thin and located at the planes z n = nd (n 
is an integer). The vector potential A due to the pan- 
cake vortices at superconducting layers satisfy a following 
equation: 



-AA- 



A 



2ttA 



^5(z-z„)^ p V( 2 Vn(r-p). (32) 



The currents in equation Li 21) are induced by pancake 
vortices with the vorticity S p = ±1 placed at the posi- 
tion p in the n-th plane. The Coulomb gauge V • A = 
was used. In addition, A z — because the direction of 
V^Vn is along the layers. It is useful to introduce an 
auxiliary potential A(r,z) = S(z — z„VA(r, z) con- 
fined to the layers, the "London vector" |l6( <f> n ( r ) = 
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J2 P < 5pfj^ 7< ' 2 Vn( r— p) and corresponding auxiliary vector 
(f> n (r,z) — ^2 n 4> n {r)5(z — z n ). In terms of these variables 
equation l)32|) can be rewritten as follows: 



AA + -A 



A 2ttA 



(33) 



Equation (|33|l can be solved by Fourier-transformation. 
The partial result is: 

Afa»)-S«-*- Mjfcffi> - (34) 

where A(q, k) is the Fourier-transform of the vector- 
potential A(r, z) and A„(q) is the plane Fourier- 
transform of the vector-potential A(r, z n ) taken at the 
n— th superconducting plane. Expressing this value by 
the inverse Fourier-transformation, we find a system of 
equations for A„(q): 

^ ( 2A^ e ~ 9 ' m ~"' d + <5 " m)A " (q) 



,— q\m— n\d 



2Aq 



(35) 



We apply eq. (|3*5j) to study the simplest case of two 
superconducting layers. Let only one pancake vortex to 
be placed in the center of the layer z = at p = 0. The 
other layer is located at z = d. The solution of the system 
(|35|l for this situation reads: 

A i(q) = i , 7tz i ,A2„ 2 — -w^i(q) 



A 2 (q) = 



1 + 4Xq + 4A 2 q 2 - e" 2 ^ 
2A qe - 2qd 



1 + iXq + 4X 2 q 2 - e- 2< i d 
In the limit qd <C 1 this solution becomes simpler 

1 



01 (q). (36) 



Ai(q) 



-0i (q) 



2 + 2A<? 

A2(q) = 2rW Mq) - 

The current density in each layer is given by: 



(37) 



Ji(q) 
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J2(q) = "4^A A2(q 



(38) 



Returning to the space coordinates, we find from eqs. 
(13 71) . (|38|) in different asymptotic regions: 



r > A d < r < A 

J n r ) le^ArV 8^A7^ 
J2(Jj 16Tr 2 Ar ^ 4A2V 



(39) 



The force between two pancake vortices isF = — *zx J, 
where J is the current produced by one pancake at the 



center of another one. It follows from eqs. I|39l ) that 
the energy of interaction between two pancakes with the 
same vorticity at the same layer is logarithmic and repul- 
sive at the distance R 3> d. A peculiarity of the few-layer 
structure is that the interaction energy of two pancake 
vortices with the same vorticity at different layers, sepa- 
rated by the distance at R ^> A, is logarithmic and at- 
tractive. It has the same absolute value as the repulsion 
of the pancakes in the same layer. This result dramat- 
ically differs from the interaction energy of two vortices 
at different layers, when the total number of layers is 
infinite: in this case the interaction in different layers 
differs from the interaction in the same layer by a small 
pre-factor d/X . The result on the logarithmic attraction 
of two pancakes in different layers and its amplitude we 
have obtained persists at any number of layers N pro- 
vided Nd <C A. It can be interpreted as the attraction 
of two "half- vortices" in the two plane, one carrying the 
flux +(f>o/2, other carrying — </>o/2. 

In the two-layer system the asymptotics for the compo- 
nents of magnetic field produced by the pancake vortex 
located in the plane z = at its origin directly follow 
from equation i|37|) . In the range r> A they are: 



B, 



M\z\ + \z~d\) 
87r(z 2 + r 2 ) 3 / 2 



(40) 



B, = J^-sgn(z)(l-^L=)f 



8-irAr 

00 



\/r 2 + z 
sgn(z — d)(l _ )f 



8nAr 

0o (2z - d) 



\/r 2 + z 2 



87r(r 2 + z 2 ) 3 / 2 ' 
In another asymptotic region d « r « A wc find: 



B, 



AkA\Jt 2 + z 2 



B< 2) (r,z) 



47rA; 



-sgn(z)(l 



Vr 2 + 



(41) 

(42) 
(43) 



Due to the strong screening effect exerted by one layer 
to another the magnetic field decays more quickly in the 
z-direction than in r -direction. The total magnetic flux 
through the plane z = and z = d is <f>(z = 0) = B vz (q = 
= °) = ^R^o « 0o/2, and $(z = d) = B vz (q = 
0, z = d) = 2 A+d 0° ~ 0o/2- The two fluxes are not 
exactly equal, and the net flux d/(2A+d) escapes through 
the remote side surface. 

The self-energy of a single pancake vortex reads: 



1 



8ttA 
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8^nV 



[|<Mq)| 2 -<M~q)-A 1 (q)] 



I6ttA £ 



d 2 q 

d 2 q ,4>l _ 4>l 1 

(2tt) 2 V 2g 2 (l + A g ) J 
R s 



(44) 
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Here R s is the lateral linear size of the sample. Due 
to divergence of E sv , it is energy unfavorable to produce 
single pancake in a layer below the Berezinsky-Kosterlitz- 
Touless transition. The energy of a pair of pancake vor- 
tices located one opposite another at different planes is: 



lv 



d 2 q 



8ttA J (2tt) 2 
•(A t ,i(q) + A„ 2 (q)] 
1 f d 2 q ,<fl 



[|^i(q)| 2 -^i(-q) 



4ttA J (2tt) 2 V q 2 (l + Aq)' 



4 ln A 

8^A £ 



(45) 



The interaction of two such a pairs separated by a dis- 
tance R ^> d is: 



A[l<Mq)(i + ^ q ' R )l 2 



8ttA J (2tt) 2 ' 
-&(-q) • (A„!(q) + A„ 2 (q)|l + e -*»- R | 2 ] 
-2Ei v 

4>l [ MqR) 



Vu(R) 



~ An 2 J 1 + Aq dq 

In the last step we used the formula |l7| : 

— ^— J (cx)dx = J[H (cz) - N (cz)} , (47) 
x + z 2 

where Hq(x) is the zeroth Struve function, and Nq(x) is 
the zeroth Neumann function. The asymptotics of the 
interaction energy l|4t)[) at small and large distances are 
as follows: 



Vu(R) 



A^i (r<A) 



i7T 2 R 



(r > A) 



(48) 



Equation (|35|) can be solved by the same method for 
any number of layers, though calculations become more 
cumbersome. However, in the region R 3> Nd equation 
(|35|l can be solved quite easily. The vector potential of a 
line vortex, identical at all layers reads: 



A x 



A 



N 



iN(f>oq x z 
q(N + 2Aq) 



(49) 



Equation H49|) allows to calculate the magnetic field , the 
current, and the interaction energy. Specifically, the sin- 
gle line self-energy and the interaction energy of two lin- 
ear vortices for an AMayer superconductor are: 



Eu, = 



A 

16tt 2 A £ 



V U (R) 



N<t>l 
4tt 2 „A 



in 2 R 



In A ( r <A) 
(r > A) 



(50) 



(51) 



From these equations we see that the energy of the vortex 
line in a few-layer system is the same as that of the Pearl 
vortex in a thin- film superconductor (if we replace Nd by 
d s ), but their interaction at short distances is N times 
stronger than the corresponding Pearl vortex interaction. 
At long distance, the interaction energy is the same as 
for the Pearl vortices. 

Next, we discuss ferromagnetic textures in a few-layer 
system. We assume that the SC and FM layers form 
very thin bi-layers separated by a finite distance d. The 
London-Pearl equations for the vector potential A m in- 
duced by the magnetic layers and screened by supercon- 
ducting layers are: 



AA, 



A 



47T V" V x [mS(z - z n )\ 



(52) 



Comparing it with equation (|32J) . we find that they be- 
come identical if we replace i4>oq x z/q by i47rm g Aq x z. 
Therefore, it is straightforward to obtain the result for 
the magnetic vector potential from that for the vec- 
tor potential induced by superconducting vortices. The 
Fourier-transform of vector potential at each layer pro- 
duced by an FM texture, identical in each plane, reads: 



iANirnigAq x z 
N + 2Aq 



(53) 



Equations l|49|) and (|53fl allow to calculate the interac- 
tion of ferromagnetic textures and vortex-ferromagnet in- 
teraction energy given the configuration of the magnetic 
texture. 

Let us consider the spontaneous stripe structure in a 
few-layer ferro-superconducting system. Under the same 
assumption about the stripe width L ^> A and the av- 
erage distance R between vortices R ^S> A, we find from 
equations (|50l) and H51|l that the interaction energy be- 
tween two vortex lines is the same as that in single layer, 
but the single line energy increases N times. The total 
vortex-ferromagnet interaction energy also increases N 
times, because the magnetic vector potential and vortex 
vector potential both increase N times if L ^> A. That 
means that the condition m4>o > e v required for spon- 
taneous formation of vortices and anti-vortices does not 
change. The domain width for a few-layer is: 



L' s ^4 ex P(^2 -C+l). 
6 4 F UNm 2 ' 



(54) 



The factor 1/N in the exponent (|54|l significantly reduces 
the domain width in the few-layer system. The sum of 
widths of parallel and antiparallel domains in an external 
magnetic field (the period of the domain structure) is: 

22/ 

L'(B ext ) = ; (55) 



1 



V 2NTvm I 
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whereas the ratio of the width parallel domain to the 
period reads: 



t = arctan — — — — 



(56) 



The critical field, at which the stripe structure vanishes 
follows from eq. 1)55(1: 



TDCl 



(57) 
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Note that it is proportional to the number of layers. The 
shift of the transition temperature AT^ is: 



AT: = 64 exp( 



am.r 



ANm 2 



It also grows with increasing TV. For the case of few SC 
films with the square array of of FM columnar dots, the 
shift of the SC transition temperature is the same as for 
a single SC film with the square array of o ferromagnetic 
dots. This result can be readily seen from the observation 
observe that A 3> Na near the transition temperature. 
Then equation l|50|) implies that the vortex line energy in 
a few-layer system is N times larger than for one layer. 
The total Ginzburg-Landau free energy for few layers is 
N times larger than that for one layer l(28|l . Thus, equa- 
tion for n s will not change and the shift of the transition 
temperature will not change either. 



6. CONCLUSIONS 

We studied the SC transition temperature in hetero- 
geneous FM-SC system in the London's approximation. 
The stripe structure of FSB, if exists, leads to the posi- 
tive shift of the transition temperature.lt can reach the 
value AT C ~ 0.1T C , when the width of the stripe becomes 
comparable with the effective penetration depth. The- 
ory predicts that the shift of the transition temperature 
is proportional to the square of magnetization per unit 
area m 2 . This fact can be checked experimentally. It was 
demonstrated that that the stripe structure must vanish 
at a very small external magnetic field between 1 and few 
tens Oersted. Simultaneously the transition temperature 
changes by the same value AT C ~ 0.1T C . This theoretical 
prediction opens way to may be the strongest magneto- 
resistive effect. 

In the multilayers the value of the value of the magnetic 
field at which the stripes disappear increases proportion- 
ally to the number of layers and the shift of the transition 
temperature grows even faster. 

The shift of transition temperature in the supercon- 
ducting layers supplied with a periodic array of magnetic 
dots may be of the same order of magnitude, but is neg- 
ative at reasonable values of parameters. 
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